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0\ '. Abstract 



We analyze the spectrum of the "local" Iwatsuka model, i.e. a two- 
dimensional charged particle interacting with a magnetic field which is 
homogeneous outside a finite strip and translationally invariant along 
'^ I it. We derive two new sufficient conditions for absolute continuity of 

the spectrum. We also show that in most cases the number of open 
O ■ spectral gaps of the model is finite. To illustrate these results we 

investigate numerically the situation when the field is zero in the strip 
being screened, e.g., by a superconducting mask. 



1 Introduction 

There has been a renewed interest recently to problems of magnetic transport 
in two-dimensional systems. Several papers ||MMP| , |FGW1| , pP|| investigated 



the influence of a weak disorder on current-carrying edge states in a halfplane 
Ha], |MS||, or in a more general domain |[FGW2||: roughly speaking, it has 



been shown that away of the Landau levels the transport survives for both 
the confining-potential and Dirichlet-boundary border of the halfplane, and 
that the same is true for Dirichlet regions containing an open wedge. 



On the other hand, it is known for long that transport can be achieved 
by a suitable variation of the magnetic field alone |p!w|| if the latter is trans- 



lationally invariant in one direction. A simple example is represented by 
the field assuming two different nonzero values in two halfplanes. It has an 



illustrative semiclassical interpretation in terms of cyclotronic radii ||CFKS 
Sec. 6.5]. On the other hand, the spectrum of this generic Iwatsuka model 
is purely absolutely continuous, while classically there are many localized 
states; this shows that the extended character of quantum states plays here 
an important role. 

Furthermore, the transport does not require the existence of different 
asymptotics; it is possible even if a variation of a constant nonzero magnetic 
field is restricted to a planar strip, i.e. its transverse profile has a compact 
support. Two sufficient conditions for absolute continuity of the spectrum 



are known in this case. In the original paper |p.w|| it was shown to be the case 



if the field derivative has opposite signs to both sides of a given point. Later 



Mantoiu and Purice ||MP|| derived what could be called the "first and third 
quadrant" condition, i.e. the requirement that ±B{xi) < ±B{xo) < ±B{x2) 
holds for some xq and all xi < xq < X2. 

In the next section we discuss the "local" Iwatsuka model. For the sake 
of completeness we describe briefly the procedure of its solution; we skip 



some details because the subject is discussed thoroughly in |w|, [MP|] . We 



shall derive a pair of new sufficient conditions for the absolute continuity. 
The first one is a global assumption of the above type: it requires that 
the field variation is nonzero and does not change sign in the strip under 
consideration. The other is a rather weak local condition: it is enough that 
the variation has at least a powerlike growth from zero at the vicinity of the 
strip edge. The main new ingredient in the second case is a semiclassical 
estimate of the transverse eigenfunction tails. Notice also that we shall not 



need the assumption of global field positivity employed in ||Iw| , |MP|| . Our 
conditions strengthen the existing results and represent a step towards the 
proof of a conjecture put forth in |PFKS| , Sec. 6.5] which states that any 



nonzero (translationally invariant) field variation spreads the Landau levels 
into a purely absolutely continuous spectrum. 

One is also interested, of course, how the spectrum looks as a set. In 
particular, it is natural to ask whether the spectral bands can overlap, and 
if so, what is the number of open gaps left. This question is addressed in 
Section 3. We show there that the number of open gaps is finite provided 



the field variation has a nonzero mean. The case of zero mean is more 
delicate and we limit ourselves to demonstrating that the conclusion persists 
under a particular assumption. We believe, however, that the restrictions are 
rather of a technical nature and we conjecture that the finiteness of open-gap 
number holds generally for the local Iwatsuka model, up to possible regularity 
requirements imposed on the field. 

The above mentioned mathematical properties of the model represent 
just one of possible motivations. On the physical side, there has been an 
interest recently [[RPM|| to the behaviour of an electron gas in a thin film 
under influence of a perpendicular magnetic field which is screened in a part 
of the plane by a superconducting "mask" . A circular "anti-dot" represents a 



solvable case which is treated in Ref. |[RPM|| within several different models 



Suppose now that the mask has the form of a straight stripQ. Then 
we get a particular case of our model, and consequently, the magnetic field 
creates states which are localized transversally and move along the strip. 
The form of these states, in particular, their "distance" from the strip is 
controlled by the velocity; the transverse eigenenergies of fast states are close 
to the Landau levels of the constant "external" field. In Section 4 we shall 
investigate numerically two variants of such a system, one with an ideal 
screening and one with an "overshoot": we will compute the band profiles 
and illustrate the band overlapping.]^ 

2 Local Iwatsuka model 

2.1 Description of the model 

We consider a charged particle confined to a plane and interacting with a 
magnetic field perpendicular to the plane. Our basic assumption is that field 
is translationally invariant in the y-direction, nonzero and constant away of 



^A negative image situation with a field constant within a strip and zero otherwise 



was treated in |Ca|. Here the spectral picture is more complicated: there are transverse 



eigenvalues of "slow" enough states embedded in t he continuum. The same system and 



several related structures were studied also in [ PM |. 

^ After finishing this paper we learned about recent results of Kim et. al. KISC j who 
studied properties of a nonmagnetic ring. In that case, of course, the spectrum remains 
to be pure point but one can see a transport analogous to the one discussed here in form 
of "edge states" moving along the ring. 



a strip of a width 2a: 

(a) the functional form of the field is 

B{x,y) = B{x) = B + b{x), 

where B > and b is bounded and piecewise continuous with supp6 = 
[—a, a]. With an abuse of notation, we employ the same symbol for functions 
on M and R^ if they are independent of one variable. 

As usual in translationally invariant situations we employ the Landau gauge, 
i.e. we choose the corresponding vector potential in the form 

A^ = 0, A.y{x) = Bx + a{x) 

with 

a{x) := / b{t)dt 



JO 

We adopt the natural system of units in which 2m = h = c = \e\ = 1 and 
assume also that the particle charge equals — |e| having in mind an electron. 
Then the Hamiltonian H = H{B, b) of our system is 

with the appropriate domain in L^(]R^). 

Since H commutes with y-translations, it allows for a standard decom- 



position. Mimicking the argument of |^, Sec. 2] we find that it is unitarily 
equivalent to the direct integral 



H{p)dp (2.1) 

with the fiber space L^(M) and fiber operator 

H{p) = -dl + {p + xB + a{x)y (2.2) 

Since the function a is bounded, the potential is for a fixed p &M. dominated 
by the oscillator term, D{H{p)) = D{—&1) fl D{x^), the spectrum of H{p) is 
purely discrete and consists of a sequence of eigenvalues e„(p) accumulating 
at +00 — see ||RS| , Thni XI11.16]. 



In the absence of the perturbation b the spectrum consists of the Landau 
levels, {(2ra + 1)B : n G Nq}. We shall first show that the latter belong to 
the spectrum in the perturbed case too, at least as its accumulation points. 



Lemma 2.1 e„(p) -^ {2n + l)B as \p\ -^ oo for any n G Nq- 
Proof: We introduce a new variable z and numbers a^ , a^ by 

z = x+p/B, a^ = b{x)dx. (2.3) 



Then we have H{p) = -dl + Vp{z) with Vp{z) = B\z + B'^ a{z - p / B)f . 
Outside [p/B — a,p/B + a] the potential V^(-2) is quadratic and equal to 
{z + a^Y for z < {p/B — a) and z > {p/B + a), respectively. Thus Vp 
converges for \p\ ^ oo pointwise to the potential of the harmonic oscillator 
Hamiltonian, 

Ho = -dl + Vo, y^{z) = Bz^ . (2.4) 

Now we employ a simple trick ||BEZ|| : we check the resolvent convergence on 



functions / = {Hq — /i)0 with (p G C^(]R). For any /i G p{H{p)) we have 

\\{H{p) - fi)-'f - {Ho - /i)-V|| = \\{H{p) - f.i)-\Vo- Vp)<p\\ ^ 

as IpI ^ oo in view of the compact support. However, such / form a dense 
set in L^(M), and therefore H{p) -^ Ho in the strong resolvent sense and the 
claim follows by Theorem VIII. 14 of [|Ka)|. | 



2.2 Absolute continuity 

Next we want to give sufficient conditions under which the spectrum is ab- 
solutely continuous. We start with analyticity of the perturbation. 

Lemma 2.2 {H{p) : p G MU {oo}} is an analytic family of type (A). In 
particular, each e„(-) is an analytic function. 

Proof: Let us first check the analyticity at infinity. We write 

H{p) = Ho + a{x){2p + 2Bx + a{x)) 



where Ho is oscillator Hamiltonian (|2.4| ). Since a is bounded by assumption, 
we have to check ifo-boundedness of the second term at the r.h.s. with a 
suitable relative bound. From the Schwarz inequality we have 

Ik/ll < ll(/ + zY'4 \\{i + ^')/ll < c{\\z'f\\ + 11/11) (2.5) 



since z h-> z{\ + z^) ^ belongs to L^(R). Given / G D{z"^) denote frjiz) : = 
rf/'^firfz)] then Hj;™"/,,!! = r]^"™!!^™'/!!. Thus (|2.5| ) acquires the form 

\\zf\\<Cr'\\z'f\\+Cr^\\f\\; (2.6) 

in other words for any a' > we can find /?' > such that 

11^/11 <«'||^Vll+/3'll/l|. (2.7) 

On the other hand, 

||(P^ + 5V)/f = {{P^ + B^z" + 2B^Pz^P + B^[P,[P,z^]])fJ) 
= llPVf + B%z^f\\^ + 2B^zPf\\^ - 25^11 ^112 



with P := —idx holds for all / G i5(]R), which yields the following inequality: 

\\z^f\\<B-^{P^ + B^z')f\\ + V2B-'\\f\\. (2.8) 

Combining (|2.7|) and ( p.8|) we arrive at 

\\zf\\<a\\{P' + B'z')f\\+(3\\f\\ (2.9) 

where a can be made arbitrarily small and [3 = [3' + a' B^^\/2. To check the 
analyticity at any point p we employ the inequalities 

\\H,f\\ < \\H{p)f\\ + \\a{-){2p + 2B-+a{-))f{-)\\ 
< ||iJ(p)/|| + 25||a|Ulk/|| + ||a"' 



loo 



Applying (p?6|) to the second term we get 

(1 - 2a'S||a||oo)||i^o/|| < \\H{p)f\\ + {2(3'B + l|a|U)l|al|c 

and since a' can be made arbitrarily small, the sought if (p)-boundedness 
of H{p') — H{j)) follows from ( |2.9| ). To conclude the proof, it is sufficient to 
apply |K|], Thms VII.2.6 and VII.1.8. | 



I.e. 



Let iljn{-,p) denote the eigenfunctions of the operator ( p.2| ) 

H{p)ll)n{x,p) = en{p)lpn{x,p). (2.10) 

Recall that without loss of generality ipn can be chosen real- valued. Following 



|[lw|| we put 

Qn,p{x) = {p + Bx + a{x)f - tnip) (2.11) 
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so ( f^.lOl ) reads ii^{x,p) = Qn,p{x)ilJn{x,p), and 

ln{x,p) = ij'„{x,pf - Qn,pi^n{x,pf. 



Lemma 2.3 (|Iw, Lemma 3.4]) ln{x,p) ^ as \x\ -^ oo. 
Finally, we have to know that the function 

fn{x,p) := (p + xB + a{x))^n{x,pf 

which determines the potential part of the energy form decays fast enough 
for large p. 

Lemma 2.4 For any p large enough there is c{p) > such that 



5c(p) e-P(^-^o) > fn{x,p) > £Me-3p(^-^o) 



f2.12) 



holds for all —a < Xq < x < a. 



Proof: By Lemma ^]T] there is po > such that Qn,p{x) > holds for all 
X G [—a, a] and \p\ > Pq. More than that, Qn,p{x) grows for a fixed x as 
\p\ —>■ oo, which makes it possible to to employ a semiclassical form for the 
tails of the eigenfunctions: for large enough p we have 



where 



and 



Mx,p) = ^il exp 

^n,p \X J 

\qn,pix)\ < exp 



XO 



Q„,p(OdeWl + g„,p(x)) 



\KA^')\dx' 



XO 



1 



(2.13) 



(2.14) 



FnA^) = J [QnTi^)^ (Qn,f (^)) | ^X 

by foil, Thm VI.2.1]. Substituting for Qn,p{x) from ( p.ll|) , we get 

5{p + xB + a{x)f{B + b{x)f 



^n,p\-^) 



4[(j9 + x5 + a(x))2-e„(p)]'/' 

1 (5 + b{x)Y + h'{x){p + xB + a{x)) 

2 [(p + xB + a{x)y ~ e„.{p)f^ 



Thus the integrand in ( 2.14 ) can be made arbitrarily small for large enough 
p. Consequently, to a fixed A > 1 we can always find px such that 



exp 



\FUx')\dx' 



< X. 



(2.15) 



holds for all p > px. The representation ( 2.13| ) is valid at the halfiine 
X > Xq, so the coefficient Ci{p) is nonzero; without loss of generality we 
may suppose that it is positive. The behaviour of ipn{x,p)/4'n{xo,p) is for 
X, xq G [—a, a] determined essentially by the exponential factor, because 
{Qn,p{xo)/Qn,p{x)) ' can be then included into the error term. Since a(-) 
is bounded and we consider x, Xq with a limited range, one has 



-p < 



Qn,piO < 2^ 



for all p larger than some pi > 0, and therefore 



2p 



6p 



if p > max(pi,p3/2), where P3/2 refers to A = 3/2 in (|2.15| ). To conclude the 
proof, notice that limp^oo(p + xB + a{x))p^^ = 1 holds for x G [—a, a] and 
put c{p) := ci(p)2. I 

Remarks 2.5 (i) For the sake of simplicity, we use large enough p in the 
above claim. However, the assumption (a) is not sensitive to mirror transfor- 
mation which changes p to —p, so the analogous result is valid as p ^ —00. 
(ii) We need both an upper and a lower bound, hence we cannot employ 



Hel as the authors of Ref. BPl did; we have to resort to a more tradi- 



tional semiclassical method. 

Now we are in position to prove the announced result under one of the fol- 
lowing additional assumptions: 

(b) b{-) is nonzero and does not change sign in [—a, a], 

(c) let ai < Gr, where we have put ai := sup{x : b{x) = in (—00, x)} and 
ttr := infjo: : b{x) = in (x, C)o)}. There exist Co,6 > and m eN such that 
one of the following conditions holds: 



\b{x)\ > Cq{x — ai)'"^ , xE[ai,ai + 6) 



(2.16) 



or 



\b{x)\ > Co{ar — x)"^ , xE{ar — 6,ar]. 



(2.17) 



Theorem 2.6 Assume (a) and (b), or (a) and (c); then \e'^{j))\ > for each 
n G No and all \p\ large enough. In particular, the spectrum of H is absolutely 
continuous. 



Proof: To prove the absolute continuity, it is sufficient by Thm XIII. 86 of |PS 



to show that £«(■) is not constant for any n G Nq. The Feynman-Hellman 
formula implies 

en(p) = [i'n{x,p),— ijJn{x,p)] =2 {p + Bx + a{x))ijjn{x, pfdx 

(2.18) 
Let us ffist investigate the integral on the semi-infinite intervals (—00, —a] 
and [a, 00). Since 



l'^ix,p) = -2{B + b{x)){p + Bx + a{x))iJn{x,pf 
and b{x) = for |x| > a, we can write 

2 / {p + Bx + a{x))^nix,p)'^ dx 

J{—oo,—a]U[a,oo) 

l'n{x,p) dx = — ln{a,p) - lni-a,p) 



(2.19) 



B 



(— 00,— a]U[a,oo) 



B 



where we have employed Lemma |2.3| . Using (|2.19|) for the second time, we 
can rewrite the last expression further as 



2 
B 



[B + b{x)){p + Bx + a{x))ilJn{x,pY dx 



and thus Eq. (|2.2|) acquires the form 



(^'nip) = ~~^ I b{x){p + Bx + a{x))^nix,p)'^dx 



■^ I b{x)fn{x,p)dx. 



(2.20) 



The first claim of the theorem follows immediately since /„(x, p) has a definite 
sign in [—a, a] for \p\ large enough. 

Assume now that (c) is valid. We shall consider the condition ( p^.l6[ ) and 
suppose that b{x) > in (a^, ai + S); the other cases can be treated in the 
same way (cf. Remark |2.5| ). In view of Lemma T^ the r.h.s. of (2.2U) can 
be estimated as follows: 



b{x)fn{x,p)dx + b{x)fn{x,p)dx 



r(n) f"-t+o par 

^^ / 6(x) e-3p(^-"^) dx - 5c(p) / |6(x)| e^^^^-"^) dx 



> 

~ 7 / 



I .la, -I ai+5 



>\c,c{p) [ re-'^«de-10acoc(p) 



'0 

Estimating the exponential function in the first integral at the r.h.s. from 
below by max{0, l — 3p^} we get 

<<^)<-^{ JTiKr.2) -^°"W--i<° (-) 

for all sufficiently large p. | 

we 



Remark 2.7 Since the perturbation is analytic at infinity by Lemma p72 
know in fact more than ( |2.21D : under the assumptions of the theorem there 
are nonzero c± and positive integers m^^^ such that 

e„(p) = {2n+l)B + c^p-""'^^ + 0{p-"'^^^~^) 
as p — > ±oo. 



3 Number of gaps 

In addition to the absolute continuity given by Theorem |2.6| we want to know 
how the spectrum of H looks like as a set. It follows from ( ^.11 ) and ( |2.2| ) 
that cr(i7) consists of a union of spectral bands /„: 

inf e„(p), sup e„(p) 



the question is how many gaps between them remain open. In this section we 
will show that their number is finite because the bands overlap at sufficiently 
high energies. 
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3.1 Field variation of a nonzero mean 

We shall distinguish two different cases depending on whether the functional 
A[b] := j_ b{x) dx vanishes or not. Suppose first that A[b] is nonzero. In 
that case we have: 

Proposition 3.1 Assume j_ b{x) dx 7^ 0. Let n{E,p) and %(£') be the 
numbers of eigenstates of H{p) and Hq, respectively, with the eigenenergy 
smaller than E. Then for any m eNq there exist pq and E{m,pQ) such that 

(rao(-E) — n{E,pQ)) sgiaA[b] > m 

holds for all E > E{m,po). 

Proof: The assumption j_^ b{x) dx 7^ is equivalent to a~ 7^ a"*". Suppose for 
definiteness that A[b] < 0, i.e. a~ > a^ . Since we are interested in the high 
energy limit we may accept without loss of generality that the field variation 
support lies in the classically allowed region, E > [po + xB + a{x)Y for any 
X G [—a, a], and to employ the Bohr-Sommerfeld quantization condition: by 
i|, Thm 7.5] we obtain then 

l-x-{E) 

7m{E,po)= ^E - {po + xB + a{x)y dx + C(E°) (3.1) 

Jx'^(E) 



-x^iE) 
'x'iE) 

where the classical turning points 



x^(E) = -^+^°-"", x^iE) = y^^^^f^^, 

satisfy by assumption the inequalities x^{E) < —a and x'^{E) > a. 



The idea is to compare (|3.1| ) with the analogous expression for Hq. Since 
the spectrum is not affected by a shift of the potential, we change the variable 
in (|2.4|) as follows, z —>■ z + a~ /B, and obtain 



7Tno{E)= ^E - (po + xB + a-y dx + 0(1) . 

JxUE) 



We have used the fact that by construction the left turning point is the same 
for both potentials, whereas the right one is moved to 
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Since a~ 7^ a"*" we have a < Xq{E) < x^{E). Taking further into account 
that the two potentials coincide to the left of —a, we may write the sought 
difference as 



7r[niE,po)~noiE)]= f \^^/E - (po + xB + aix)^ 

-^E- {po + xB + a~)^j dx+ ^^E-{po + xB + a+y 

, T /•-'■(^) , 

- y/E - {po + xB + a-y\dx+ / ^yE - (po + xB + a+)2 dx 

^ JxUE) 



xl{E) 

+0(^°), (3.2) 

the last term being simply a positive number independent of E. 

In the first term at the r.h.s. of ( |3.2| ) we integrate over a fixed interval, 
hence the result is 0{E~^/'^) as E ^ 00 and may be absorbed into the error 
term. Furthermore, choosing po ^ ~o,^ ~ o.B we achieve that the integrand 
in the second term is positive, hence we have 



7r[n(E,po)-no(^)] > / ^E - (po + xB + a+f dx + C>(£;°) . 

JxUE) 



It remains to estimate the last integral. Since the function is nonnegative, 
decreasing and vanishes only if x = x^{E), we take any 5 G (0, x^{E)—Xq{E)) 
and use a simple bound. 



px''{E) i'x''{E)-5 

/ y^E- {po + xB + a+ydx > / ^/E- {po + xB + a+y dx 

JxlUE) JxUE) 



a — a~^ 



> \/2B6VE-B^6'^l 5 

which yields the sought result for E large enough. The inequality 

n{E,po) < no{E) - m 

for A[b] > is obtained in the same way. | 

Corollary 3.2 If A[b] 7^ the number of open gaps in the spectrum of H is 
finite. 
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Proof: Let again a > a~^. Since e„(p) — > (2n + 1)5 as p — > oo for a fixed 



n G N by Lemma |2.1| , it is sufficient to find n and p such that 

en+i{p)<{2n + l)B (3.3) 



holds for all n > n. This follows immediately from Proposition |3.1| with 
m = 2. In the opposite case, a~ < a"*", the inequality ( |3.3|) is replaced by 
e,_i(p)>(2n + l)5. I 

3.2 The case of zero mean: an example 

If A[b] =0 the situation is more complicated since the two potentials differ 
only in a subset of the interval (—a, a). We restrict ourselves to an example. 
As above, put ag := supja; : b{x) = in (— oo,a;)}, and suppose that there is 
a number c G (a^, a) such that a{x) < holds in {ai, c) and a{c) = . 

Let us show that the above conclusion about the finite number of gaps 
persists in this case. We employ again the Bohr-Sommerfeld condition ( |3.1|) 
choosing E and p in such a way that c will be the right turning point, 
E = {p + Bc)^. It may happen, of course, that there is another classically 
allowed region to the right of c but changing the potential to E there certainly 
does not increase the number of bound states, i.e. 

Tm{E,p)> ^/E-{p + xB + a{x)y dx + 0{E^) . 

This has to be compared with the number of oscillator states corresponding 
to a = 0. Introducing the variable y := c — x we get 

rc-at I 

7i[n{E,p)-no{E,p)] > / a/2/E - By + a{c-y) ^ By - a{c-y) dy 

Jo 

- / \j2\fE-By^Wydy + 0{E^). 

Jo 

The difference of the integrands can be written as 



4/ 



AE 



'l-^^^^^v/lF-^-Jl-llv^ 



^ (^^By-a{c-y) - /B^) + 0{E 
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The error term can be absorbed into that of ( ^.2| ), hence we get 

[^By-a{c-y) - v^j dy + 0{E^) , 

which gives the desired result since the integral is positive by assumption. 

The result can be modified for the other end of supp b, and in a similar 
vein one can treat cases where a(x) is locally positive close to an endpoint. 
We will not pursue the line further since we believe that a deeper analysis 
of the behaviour around the classical turning points would be useful, and at 
the same time, that the property in question is valid generally: 

Conjecture 3.3 The number of open gaps in the local Iwatsuka model is 
finite for any nonzero b satisfying minimal regularity assumptions such as 
those formulated in Section. 2.1. 

4 Example: a screened strip 

In this section we shall discuss in more detail the example mentioned in the 
introduction. Apart of the physical interest mentioned there, the system 
represents a solvable case which allows us to illustrate the results of the 
previous sections. 

4.1 The ideal situation 

We shall consider first the simplest situation (without "overshoot" in the 
terminology of |[RPM|| ) when the distance between the film supporting the 
electrons and the superconducting strip is negligible and the magnetic field 
is perfectly screened for |x| < a. Consequently, the field perturbation b{x) is 
of the form 

b{x) = -BQ{\x\-a) (4.1) 

with the vector potential 

A^ = Q, Ay{x) = B{x - a)Q{x - a) + B{x + a)Q{-x - a) . (4.2) 

The fiber operator in the decomposition ( p.l| ) then looks as follows 

H{p) = -dl + [p + B{x - a)e{x -a) + B{x + a)e(-x - a)f . (4.3) 
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The real line decomposes into a union of several parts in which the corre- 
sponding Schrodinger equation can be solved. Inside the strip the potential 
term in ( [4 .31 ) is constant and equal to p^, so the eigenfunctions ^„(x,p) of 
H{p) are there of the form 

ipnix, p) = A exp [unip) x]+ B exp [-Kn{p) x], X e [-a, a] , 

where K„(p) := a/p^ — £n{p)- On the other hand, for x ^ [—a, a] the function 
a is constant and equal to its boundary values; in view of ( p.3| ) they are 

a~ = —a^ = aB . 

In other words, outside the strip the potential has a parabolic shape and its 
two branches are shifted mutually by 2a. Thus we use the substitutions 

X + a+ p/B ... X < -a . . 

X — a + p/B ... X > a 

denote r = Bz^ and look for solutions of Eq.( |2.10| ) in the form 

i)n{z, p) = exp(-r/2) M„(r, p) (4.5) 

In a usual way we check that Un{r,p) satisfies the confluent hypergeometric 
equation | |A!J| |, 

ru'^{r,p)+ ( 2 -'^j "n(^'P) -aunir,p) =0 

with a := '^^ ■ Consequently, the general solution in the "parabolic 
regions" is of the form 

Mx,P) = e-''/' jCi M (a, i, r") + C2 f/ (a, ^, r 

We should not forget, however, that r in this formula comes from two different 
substitutions ( |4.4]) , and moreover, that the map 2; 1— >■ r is not a bijection. 
With the exception of the case p = the real line therefore decomposes into 
four intervals. In the outer regions the requirement of L^ integrability leaves 
the term with the U function only, while in the middle "parabolic regions" 
both functions are generally contained in the linear combination. Denote 

p 
xo = -— -asgn(p). 
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For p < this point where the potential reaches zero is to the right of the 
strip, a < xq, and we have 



ipn{x,p) 



A exp [nn{p) x] + B exp [— /t„(p) x] 
e-'^/2 ^c^ M{a, p, r) + C2 U{a, p, r)} 



^ D2e-'-/'U{a,f3,r) 



-00 < X < —a 
—a < X < a 

a < X < Xq 

xo < X < cxo 



(4.6) 



and a similar expression for p > when xq < —a. In the case when p = we 
set Xq = —a by definition, then Ci = and C2 = Di. Matching the solutions 
( |4.6| ) smoothly at the points —a, a, Xq we arrive at a homogeneous system of 
six linear equations for the coefficients A,B,Ci,C2,Di, and D2, which yields 
an implicit equation for the energy levels e„(p). 



4.2 The model with an overshoot 

The above model is, of course, idealized because in reality the magnetic field 
is never fully screened and the vector potential changes sharply near the 
edges of the screening strip. Another simplified description takes this effect 
into account by putting 



x\ 



b{x) = -BQ{a- \x\) + B\x\ 

which corresponds the vector potential 

A{x) =xBQ{\x\ -a). 

Recall that a numerical analysis performed in ||KPM|| for a screening 
suggests that a realistic field profile should lie between these two 
The fiber operator in ( |2.1| ) now reads 

H{p) = -dl + [p + BxQ{\x\ - d)f. 

Inside the strip we get the same solution as in the model without 
shoot. However, the boundary values of the function a(x) given by 
now equal to each other, 

"±770 



(4.7) 



by a disk 
extrema. 

(4.8) 

an over- 
lO) are 



lim 



h{x) dx = , 
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hence for x ^ [—a, a] we use the same substitutions to the left and to the 
right of the strip, 



z' = X + p/B and r' = B{z' 



/n2 



In analogy with Sec. O we find the overshoot solutions ipn{x,p) of ( [4 .51) in 
the form 

' D[ e-"'/2f/(a, /?, r') -oo < x < -a 

A' exp [unip) x] + B' exp [—Hn{p) x] —a < x < a 

e-"'/2 {C[ M{a, (3, r') + C'^ U{a, (3,r')} a < x < x'^ 



i)n{x,p) = < 



^ D',e-^'/'U{a,(3,r') 



Xq < X < CXD 



(4.9) 

where x'q = —p/B. Matching the functions ( [4. 91) smoothly at the points 
—a, a, Xo we arrive again at an implicit equation for the energy levels in{p)- 



4.3 Numerical results 



Using the Ansatze ( |4.6| ) and ( [4. 9] ) for the two described models, we have 
solved the matching conditions numerically for several values of the param- 
eters. The obtained band functions, denoted e„(p) and e„(p) for the case 
without and with the overshoot, respectively, are plotted in Fig. 0. 

We see that in the ideal situation (without an overshoot) the band func- 
tions e„(p) have the only stationary point at p = 0. This is true for any n as 
it follows from Eq. (|2.2(]| ) which now acquires the form 



en(p) = 2p / ^nix^pfdx. 



(4.10) 



The situation is completely different in the presence of an overshoot. Here 
the fiber operator ([4.8|) can be rewritten as a sum of Hq with the pertur- 
bation supported in [—a, a]. Consequently, the functions e„(p) exhibit in 
general many oscillations corresponding to the harmonic oscillator solutions. 
Following this argument we would expect that in higher levels these oscilla- 
tions will be partially "ironed" out for sufficiently large a. As an illustration, 
see the curves for ra = 6 in Figs. |l| and ^ 
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Figure 1: The band functions in the units of B for B = 10^, a = 0.01 
(dotted hne), a = 0.02 (dashed hne) and a = 0.03 (full line). (Left) the ideal 
situation. (Right) the model with an overshoot. The thick lines show the 
Landau levels. 



In Section |3.1| we have shown that if the functional A[b] is nonzero the 
number of open spectral gaps is finite. Since the field variation ( [4.1| ) of the 
ideal model gives A[b] = —2Ba, it follows from Proposition p.lj that for a large 
enough band index n, depending on a and B, the width of the spectral band 
In exceeds 2B, the distance between Landau levels. Moreover, applying the 
Bohr-Sommerfeld quantization rule we can easily derive an explicit condition 




Figure 2: Overlap of the spectral bands. (Left) the ideal situation: a = 0.015 
(dashed line), a = 0.02 (full line). (Right) the model with an overshoot: 
a = 0.03 (dashed line), a = 0.045 (full line). The field is always B = lOl 

under which this happens: 



'^overlap(o, -S) ~ - + 



TC 



Sa'^B 



(4.11) 



This is shown in Fig |^. In the case of the overshoot model with the field 
variation given by ( |4.7| ) we have no simple condition analogous to ( [4.111 ). 
Nevertheless, the assumptions of the example of Section |3.2| are satisfied 
here so we know that at sufficiently high energies the spectral bands overlap 
in this case too. 

On a heuristic level, this can be understood from the behavior of the har- 
monic oscillator solutions. It is well known that for large n they are strongly 
oscillating but the squared modulus smoothed over a small interval tends as 
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Figure 3: The squared eigenfunctions for n = 0, 6, and different values of p. 
The dashed hues indicate the position of the nonmagnetic strip with a = 0.03, 
the field is B = 10^ 



n ^ oo to the probability density of finding the classical oscillator at a given 
point [ BEH , Section 8.3]. The latter diverges at the turning points. Tak- 
ing for the difference between en{p) and the corresponding Landau level the 
first-order perturbation theory expression, we expect therefore this quantity 
to reach its maximum for the values of p corresponding to the vicinity of the 
turning points. Fig ^ illustrates that it is indeed the case. 

Once we have solved the matching conditions for the band functions we 
can also find the coefficients in ( [4.6|) and ( [4.9|) to obtain complete eigen- 
functions ipn{x,p) and ipnix,p), respectively. Fig |^ shows the corresponding 
probability densities for ground state and the excited state with n = 6. In 
the absence of an overshoot it can be seen that as the screened strip [—a, a] 
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crosses the center of the distribution , each particular peak of the density 
becomes broader and more smeared. On the other hand, for large \p\ when 
the "bulk" of the \ipn{x,p)\'^ support lies outside the strip, its shape is close 
to that of the harmonic-oscillator probability density. In the case of an over- 
shoot the picture is similar but the plot shows a significant deformation for 
the value of \p\ corresponding to the classical turning points. 
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